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ABSTRACT 


This  paper  shows  that  much  new  information  about  the  dynamics  of  combat  between 
two  homogeneous  forces  modelled  by  Lanchester-type  equations  of  modern  warfare  (also 
frequently  referred  to  as  "square-law"  attrition  equations)  with  temporal  variations 
in  fire  effectivenesses  (as  expressed  by  the  Lanchester  attrition-rate  coefficients) 
may  be  obtained  by  considering  Liouville's  normal  form  for  the  X and  Y force-level 
equations.  It  is  shown  that  the  relative  fire  effectiveness  of  the  two  combatants 
and  the  intensity  of  combat  are  two  key  parameters  determining  the  course  of  such 
Lanchester-type  combat.  New  victory-prediction  conditions  that  allow  one  to  forecast 
the  battle's  outcome  without  explicitly  solving  the  deterministic  combat  equations 
and  computing  force-level  trajectories  are  developed  for  fixed-force-ratio-breakpoint 
battles  by  considering  Liouville's  normal  form.  These  general  results  are  applied  to 
two  special  cases  of  combat  modelled  with  general  power  attrition-rate  coefficients. 

A refinement  of  a previously  known  victory-prediction  condition  is  given.  Temporal 
variations  in  relative  fire  effectiveness  play  a central  role  in  these  victory-predic- 
tion results.  Liouville's  normal  form  is  also  shown  to  yield  an  approximation  to 
the  force-level  trajectories  in  terms  of  elementary  functions. 


1.  Introduction. 

Even  though  combat  between  two  military  forces  is  a complex  random  process  (see 

Note  1 of  TAYLOR  and  BROWN  as  a consequence  of  the  pioneering  1914  work  of 

F.  W.  LANCHESTER^2^  (see  Note  1)  military  operations  analysts  since  about  the  end  of 

World  War  II  have  used  simplified  deterministic  differential  equation  models  to  develop 

insights  into  the  dynamics  of  combat  (see  Note  2).  Today  Lanchester-type  models  of 

quite  complex  military  systems  have  been  developed  in  the  United  States  (see,  for  example, 
(91 

BONDER  and  HONIG  ) and  require  a digital  computer  for  their  implementation.  A simple 
combat  model,  however,  may  yield  an  understanding  of  important  relations  that  are 
difficult  to  perceive  in  a more  complex  model,  and  such  insights  may  provide  guidance 

r q 1 

for  higher  resolution  computerized  investigations  (see , for  example,  BONDER  and  FARREL1  1 
and  WEISS  In  this  paper  we  will  examine  such  an  idealized  Lanchester-type  model 

in  order  to  obtain  some  insights  (specifically,  the  tradeoff  between  quality  and  quantity 


of  weapon  systems)  into  the  dynamics  of  combat  between  two  homogeneous  forces  with 
temporal  variations  in  weapon  system  effectivenesses. 


In  this  paper  we  develop  new  victory-prediction  conditions  that  sometimes  allow 
us  to  forecast  the  battle's  outcome  without  explicitly  solving  the  combat  equations 


(see  Note  3)  and  computing  force-level  trajectories  (see  Note  4).  We  obtain  these 
results  for  variable-coefficient  Lanchester-type  equations  of  modern  warfare  by  consid- 
ering Liouville's  normal  form  of  the  X and  Y force-level  equations  (see  p.  270  of 
INCE^2^  and  KAMKE^2^)  and  using  techniques  recently  applied  to  Lanchester  combat 
theory  by  TAYLOR  and  PARRY^^.  These  results  complement  and  extend  those  of  Taylor 
and  Parry^^,  and  they  show  that  the  key  parameters  affecting  a battle's  outcome,  at 
least  for  fixed-force-ratio-breakpoint  battles  with  the  initial  force  ratio  held  con- 


stant, are  the  relative  effectiveness  of  the  weapon  systems  and  the  intensity  of  combat. 

Such  results  are  not  only  important  in  their  own  right  but  also  useful  in  the  quanti- 

[44-461 

tative  analysis  of  tactics  (see , for  example,  TAYLOR  J). 
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This  paper  is  organized  in  the  following  fashion.  First,  we  review  F.  W. 
Lanchester's  classic  mathematical  model  of  combat  between  two  homogeneous  forces  and 
its  extension  to  cases  of  time-varying  fire  effectivenesses.  Next,  we  transform  first 
the  independent  variable  (time,  t)  and  then  the  dependent  variable  (force  level)  to 
obtain  Liouville's  normal  form  for  the  X and  Y force-level  equations.  Then  we 
develop  some  new  victory-prediction  conditions  from  Liouville’s  normal  form  and  apply 
these  general  results  to  two  special  cases  of  power  attrition-rate  coefficients. 
Finally,  we  discuss  the  significance  of  our  developments. 

2.  Lanchester's  Classic  Formulation. 

r 26 1 

F.  W.  Lanchester1  J hypothesized  (see  Note  5)  in  1914  (see  Note  6)  that  combat 
between  two  military  forces  "under  modern  conditions"  could  be  modelled  by  (see  Note  7) 


( dx/dt  = 

-ay 

with 

x(t=0)  = xQ, 

(1) 

) dy/dt  = ■ 

-bx 

with 

y (t-0)  = yQ, 

where  t = 0 denotes  the 

time 

at  which  the  battle  begins. 

x(t)  and  y(t) 

denote 

the  numbers  of  X and  Y at  time  t,  and  a and  b are  nonnegative  constants  that 
are  today  called  Lanchester  attrition-rate  coefficients  and  represent  each  side's  fire 
effectiveness.  Lanchester  considered  this  simple  model  in  order  to  provide  insights 
into  the  dynamics  of  combat  under  "modern  conditions"  and  justify  the  principle  of 
concentration  (see  Note  8).  We  will  accordingly  refer  to  (1)  as  Lanchester's  equations 
of  modern  warfare.  Various  sets  of  physical  circumstances  have  been  hypothesized  to 
yield  them:  for  example,  (A)  both  sides  use  aimed  fire  and  target  acquisition  times 

[54] 

are  constant  (see  WEISS  ),  or  (B)  both  sides  use  area  fire  and  a constant  density 
defense  (see  BRACKNEY ^ ^ ^ ^ ) . Other  forms  of  Lanchester-type  equations  have  appeared  in 
the  literature,  but  we  will  not  consider  these  here  (see  DOLANSKY^"^  and  TAYLOR^'*’  ^1). 

From  (1)  Lanchester  deduced  his  famous  square  law 

b{x2  - x2 (t) } * a{y2-y2(t)},  (2) 


2 


f 


which  has  the  important  implication  that  a side  can  significantly  reduce  its  casualties 
by  initially  committing  more  forces  to  battle.  It  follows  from  (2)  that 

X will  be  annihilated  « — < . (3) 

y0  7 b 

Unfortunately,  no  simple  relationship  similar  to  (2)  holds  in  general  for  variable 
attrition-rate  coefficients  so  we  consider  other  means  for  developing  (3).  As  is  well 
known,  the  X force-level,  x(t),  is  given  by 

x(t)  = xn  cosh(/ab  t)  - y~  /jj-sinh(/ab  t)  . (4) 

U U*  D 


We  may  also  deduce  (3)  by  writing  (4)  as 

x(t)  = ^{(x0-y0>/^)exp(>^bt)  + (x04y0|/^)exp(-/abt)}  , (5) 

and  observing  that  x(t)  can  become  zero  if  and  only  if  the  coefficient  of  the  increas- 
ing exponential  is  negative.  We  observe  from  (5)  that  annihilation  occurs  in  finite 
time. 

[53] 

As  H.  K.  WEISS  has  emphasized,  engagements  that  continue  until  one  side  is 
wiped  out  are  rare.  Thus,  we  see  that  a model  of  battle  termination  is  required. 
Although  we  are  well  aware  that  battle  termination  is  a complex  random  process  for 
which  it  is  by  no  means  certain  that  force  levels  are  the  only  significant  variables 
(see  Note  4),  we  assume  that  combat  ends  when  either  of  two  given  "breakpoint"  force 
ratios  is  reached.  Introducing  the  force  ratio  u = x/y,  we  have  that  these  "break- 
point" force  ratios,  denoted  as  u^  when  X wins  and  u^  when  Y wins,  satisfy 
0 < u|  < Ug  = u(t=0)  < u^  £ +“>.  Corresponding  to  a fight  until  the  annihilation  of 
one  side  or  the  other  is  the  case  in  which  u^  = 0 and  u^  = +■». 

Let  us  now  consider  such  a fixed-force-ratio-breakpoint  battle.  Introducing  the 
force  ratio,  u * x/y,  we  see  that  it  satisfies  the  Riccati  equation 


du/dt  = bu2  - a with  u(t=0)  = u^  = xq/Yq*  (6) 

Let  u+  " /a/b  denote  the  positive  root  of  the  quadratic  equation  bu2  - a = 0 and 
observe  that  du/dt  <0  *»  u < u+  (see  Figure  2 of  Taylor  and  Parry In 
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particular,  u^  < u+  ■»  du/dt(t)  * du/dt(t=0)  < 0,  whence  follows 

THEOREM  1:  Consider  Lanchester's  equations  of  modern  warfare  with  constant 
attrition-rate  coefficients  (1).  Then  X will  lose  a fixed-force-ratio 
breakpoint  battle  in  finite  time  if  and  only  if  Xg/yg  < • 

We  observe  that  (3)  is  a special  case  of  Theorem  1 corresponding  to  u^  = 0 and 

u£  » +°>.  One  result  of  this  paper  is  to  generalize  [in  a way  different  from  that  given 

by  Taylor  and  Parry (see  Section  6)]  Theorem  1 to  cases  of  variable  attrition-rate 

coefficients. 


3.  Variable  Attrition-Rate  Coefficients. 

[3  5 81 

The  pioneering  work  of  S.  BONDER1  ’ ’ on  methodology  for  the  evaluation  of 

military  systems  (especially  mobile  systems  such  as  tanks,  mechanized  infantry  combat 
vehicles,  etc.)  has  generated  interest  in  variable-coefficient  Lanchester-type  equations 
and  has  led  to  improved  operations  research  techniques  for  the  prediction  of  these 
coefficients  (see  Note  10) . Let  us  therefore  consider 


dx/dt  - -a(t)y 
dy/dt  = -b(t)x 


with  x(t=0)  = Xq, 


with  y (t=0)  = yQ, 


(7) 


where  a(t)  and  b(t)  denote  time-dependent  Lanchester  attrition-rate  coefficients. 

These  coefficients  depend  on  such  variables  as  force  separation,  tactical  posture  of 
targets,  rate  of  targer  acquisition,  firing  doctrine,  firing  rate,  etc.  (see  reference  8). 
We  will  also  refer  to  (7)  as  the  equations  for  a square-law  attrition  process,  since 
an  "instantaneous"  square  law  holds  even  when  a(t)/b(t)  is  not  constant  (see  Taylor 
and  Parry  also  references  45  and  47). 

A large  class  of  tactical  situations  of  interest  can  be  modelled  with  the  follow- 
ing general  power  attrition-rate  coefficients  (see  reference  8) 


a(t)  - ka(t-KI)lJ, 


and 


b(t)  = kb(t4C+A)v, 


(8) 
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where  A,C  i 0.  We  will  call  A the  offset  parameter,  since  it  allows  us  to  model 
(with  ii, \>  i 0)  battles  between  weapon  systems  with  different  effective  ranges  (see 
Note  11).  We  will  call  C the  starting  parameter,  since  it  allows  us  to  model  (again, 
with  y,v  5fc  0)  battles  that  begin  within  the  maximum  effective  ranges  of  the  two  systems. 
Restrictions  that  must  be  placed  on  y and  v,  which  are  not  necessarily  integers, 
are  discussed  below. 

The  above  nomenclature  is  motivated  and  possible  applications  of  our  results  are 

[3  5] 

indicated  by  considering  S.  Bonder's  ’ model  of  the  constant-speed  attack  on  a 
static  defensive  position 

dx/dt  * -a(r)y  = -aQ(l-r/Ra),Jy , dy/dt  = -3(r)x  = -Bg(l-r /Rg) Vx , (9) 

where  r denotes  the  range  between  opposing  forces,  y,v  £ 0,  and  Rq  denotes  the 
maximum  effective  range  of  the  Y weapon  system.  Range  is  related  to  time  by 
r(t)  = Rg-vt,  where  Rg  denotes  the  opening  range  of  battle  and  v > 0 denotes  the 
constant  attack  speed.  Then  the  offset  and  starting  parameters  are  given  by 

A = (Rg-Ra)/v,  and  C = (Ra-RQ)/v.  (10) 

We  observe  that  A,C  £ 0 if  and  only  if  R^  i i Rg.  By  considering  (10)  and 

Figure  1,  the  reader  should  have  no  trouble  in  understanding  our  terminology  for  A 
and  C.  In  the  model  (9)  y,  for  example,  is  used  to  model  the  range  dependence  of 
Y's  attrition-rate  coefficient  (see  Figure  2). 

From  (7)  we  obtain  the  X force-level  equation 

d?-(^na(t)l^“-a(t)b(t)x  = 0,  (11) 

with  initial  conditions 

x(t-0)  = xQ  and  { [ l/a(t) ]dx/dt }t=Q  = -yQ, 

X Y X 

where  tg  ■ max(tg,tg),  and  tg  denotes  the  right  most  finite  singularity  of  the 

r 231 

X force-level  equation  (see  Ince  ).  The  coefficients  a(t)  and  b(t'  are  then 
positive  continuous  functions  Vt  > tg.  For  the  coefficient  (8),  we  observe  that 
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modelling  constant-speed  attack.  [Notes:  1.  The  maximum  effective  ranges  of  the  two  weapon 
systems  are  denoted  as  Rqj  and  Rg . 2.  The  opening  range  of  battle  is  denoted  as  Rq  and 
(as  shown)  Rq  < minimum  (Ra,Rg).  3.  The  offset  parameter  is  given  by  A = (Rg-Ra)/v.  4.  The 
starting  parameter  is  given  by  C = (Rjj-RqJ/v.I 


range  of  the  weapon  system  and  kill  capability  at  zero  range  held  constant.  [Notes:  1.  The 
maximum  effective  range  of  the  system  is  denoted  as  Ra  = 2000  meters.  2.  a(r=0)  = ag  = 0.6  X 

casualties/(unit  time  x number  of  Y units)  denotes  the  Y force  weapon  system  kill  rate  at 
zero  force  separation  (range).  3.  The  opening  range  of  battle  is  denoted  as  Rq  = 1250  meters 
and  (as  shown)  Ro  < R«.l 


tig  “ -C.  Moreover,  to  insure  the  existence  of  all  derivatives  required  in  subsequent 
analysis,  we  assume  that  the  second  derivatives  of  a(t)  and  b(t)  exist  Vt  > tg. 
We  similarly  have  the  Y force-level  equation 


a(t)b(t)y  = 0, 


(12) 


with  initial  conditions 


y(t-0)  = yQ  and  { [l/b(t)]dy/dt)t=0  = -xQ. 

It  is  necessary  to  place  further  restrictions  on  a(t)  and  b(t)  in  order  to 
insure  that,  for  example,  the  transformation  introduced  in  the  next  section  is  well 
defined.  Thus,  we  assume  that  the  following  condition  holds. 


t t 

CONDITION  (A):  | a(s)ds  and  j b(s)ds  are  bounded  for  all  finite  t tg. 

fc0  fc0 
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This  condition  also  guarantees  (see  Theorem  6. A. 2 on  p.  226  of  HILLE  ) that  (7)  has 
a continuous  solution  for  all  finite  t i 0 i tg  (see  Note  12) . If  Condition  (A)  is 
to  hold,  then  for  the  general  power  attrition-rate  coefficients  we  must  have  u,v  > -1. 


4.  Transformation  of  the  Battle's  Time  Scale. 
Let 


/a(s)b(s) ds  , . 


(13) 


and  denote  r(t“0)  as  tg.  It  follows  that  ig  i 0 for  tg  * 0.  Condition  (A) 

guarantees  that  r ■ x(t)  is  well  defined  (i.e.  bounded  for  all  finite  t £ tg)  by 

f 21 

the  Cauchy-Schwarz  inequality  for  integrals  (see  p.  123  of  BELLMAN  J)*  The  transfor- 
mation is  invertible,  since  dt/dt  ■ /a(t)b(t)  > 0 for  t > tg.  Thus,  we  may  consider 
that  t “ t(x).  We  assume  (and  give  below  conditions  that  guarantee)  that 
lim  i(t)  ■ +“  so  that  the  range  of  x(t)  is  [0,+°°)  for  te  (tn,+“). 

Considering  the  constant  coefficient  result  (4),  we  will  call  the  quantity 
/a(t)b(t)  the  "intensity  of  combat"  (see  also  Taylor  and  Parry  ^^ ) ; since  the  larger 


8 


r 


it  is, the  more  quickly  the  battle  is  moving  towards  termination.  Then  t - tq  is 

related  to  the  average  intensity  of  combat  by 

t 

T ~ Ty  ' I ' ! 

0 


l 

{(l/t)  | /a(s)b(s)  dslt  ” / a(t)b(t)  t. 


(14) 


Applying  the  transformation  (13)  to  (11)  and  (12),  we  obtain 
d2x 


dr 


. r 1 d „ rb(t)l,dx 

r+{2  57 *n  £u)J}^~x  ' °* 


(15) 


with  initial  conditions 


X(T-T0) 


and 


Xq  and  { [b(t)/a(t) ]^^2dx/dT)T_T 


d^  + {2  dV*n  [bS{|}d?"y  = °* 


'o’ 


(16) 


with  initial  conditions 


y(T-T0) 


and  {[a(t)/b(t)]1/2dy/dT>T=T^  = -xQ. 


.[50] 


Taylor  and  Brownl'’WJ  have  shown  that  the  X force-level  equation  (11)  may  be 
transformed  into  a linear  second  order  differential  equation  with  constant  coeffi- 
cients if  and  only  if 

1 d 


/a(t)b(t) 


dt 


in 


b(t)_ 


CONSTANT , 


(17) 


with  the  desired  transformation  being  given  by  t - K j /a(s)b(s) ds , where  j ...ds 
denotes  an  indefinite  integral  and  K an  arbitrary  constant.  Hence,  equation  (11) 
may  be  transformed  into  such  a constant  coefficient  equation  if  and  only  if  (15)  is  a 
constant  coefficient  equation.  We  observe  that  (15)  has  a solution  in  terms  of  ele- 
mentary transcendental  functions  in  the  quasi-autonomous  case  in  w’-.ich  the  ratio  of 
attrition-rate  coefficients  is  constant  (see  Note  13),  i.e. 

b(t)/a(t)  - CONSTANT.  (18) 


■ 


Except  when  (17)  holds,  the  solution  to  (11)  Is  complex,  and  the  qualitative  behavior 
of  force-level  trajectories  has  been  difficult  to  establish.  In  particular,  one  is 
interested  in  answering  such  questions  as 

(Ql)  Who  will  "win"?  Be  annihilated? 

(Q2)  How  do  force  levels  decrease  over  time  and  how  many  survivors  will  the 
winner  have? 

(Q3)  How  do  changes  in  the  initial  force  levels  and/or  weapon  system  parameters 
affect  the  outcome?  Is  concentration  of  forces  a good  tactic? 

(Q4)  How  long  will  the  battle  last? 

We  will  now  show  how  to  answer  question  (Ql)  without  explicitly  solving  the  Lanchester- 
type  equations  (7). 

Equation  (15)  is  highly  significant  because  it  clearly  shows  that  the  course 
of  combat  depends  on  just  the  two  weapon  system  parameters:  (1)  R(t)  = a(t)/b(t), 
the  relative  effectiveness  (Y  to  X)  of  the  two  weapon  systems,  and  (2)  I(t)  = 
/a(t)b(t) , the  intensity  of  combat  (through  equation  (13),  which  relates  I(t)  to 
t).  Both  these  parameters  may  vary  over  time,  and  equation  (15)  tells  us  that  the 
nature  of  such  temporal  variations  in  relative  effectiveness  has  a significant  effect 
upon  the  course  of  combat.  Moreover,  this  relationship  may  be  more  explicitly  seen 
by  transforming  (15)  to  Liouville's  normal  form. 


5.  Reduction  to  Liouville's  Normal  Form. 


Let  us  assume  that  Tq  > 0 (i.e.  0 > tg)  so  that  a(t)/b(t) 

entiable  and  satisfies  0 < a(t)/b(t)  < +°°  for  0 * t < +°°.  Let  a ^ 

etc.,  and  recall  that  t * t(x).  The  substitution 

—1 1/ 4 

late; /a 

x (t ) * X(t) 


b(t)/b 


is  twice  differ- 
denote  a(tmO) , 


(19) 


transforms  (15)  into  the  so-called  normal  form  (see  RAINVILLE^^)  with  the  first 
derivative  of  the  dependent  variable  removed 


(1  + F(t)}X 


0, 


(20) 
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I 


with  initial  conditions 


XCt-Tq)  - x0,  and  dX/dx(x-xQ)  = -y0/a0/b0  " x0e0* 


where 


F(t)  - P"(t)/P(t), 
R(t)  - a(t)/b(t) , 


P(t)  - [R(t)f1/4, 


e(t) 


4/a(t)b(t) 


(21) 

(22) 


Eq  denotes  e(t“0),  and  P'(t)  denotes  dP/dx.  Equation  (20)  is  Liouville's 


celebrated  normal  form  (see  p.  23  of  LIOUVILLE1  J)*  Similarly, 

T 1/4 

idvc;/d 

Y(t)  = Y(x) 


yields 


with  initial  conditions 


dzY 

dr7 


a(t)/a. 


- { 1 + G(t)}Y  = 0, 


(23) 


(24) 


where 


Y(t-t0)  - yQ,  and  dY/dT(r=T0)  = -x0/b0/a0  + y0c0. 


1 / 4 

G(t)  = Q"(t)/Q(t)  , and  Q(t)  = [R(t)]A^  = 1/P(x). 


(25) 


We  observe  that 


F(x)  + G(x)  = 2(^inP)2. 


(26) 


It  is  sometimes  convenient  to  express  F and  G in  terms  of  the  old  time 
variable  t.  Then 

1 J J 1 J J JD 

(27) 


F(t)  - 7^TT^t-^^b(t)-i^rinR(t)+  , 


4bz (t)  dt  1 dt  4 dt 


and 


Observing  that  (21)  may  also  be  written  as 


F(x)  » j^yin  P(x) +{^ln  P(x)}2  , 


(28) 


(29) 
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we  see  that  d An  P/dT  - constant  (or,  equivalently,  e(t)  - constant)  implies  that 
F(t)  is  constant,  although  the  converse  may  not  be  true. 


0.0 

Writing  (20)  as  d X/dt  - X » F(t)x,  we  may  use  variation  of  parameters  (see 

r 23 1 

pp.  122-123  of  Ince  ) to  obtain  the  solution  to  (20)  as 

X(t)  - xQ  cosh  (t-tq)  - [y0/a0/b0  + x0e0]sinh(T-x0) + j F(o)sinh(T-o)X(o)do . (30) 

To 

In  terms  of  the  original  time  variable  t and  dependent  variable  x,  we  have 


x(t)  - 


a(t)/ar 


“1 1/ A 


b(t) /b. 


(xq  cosh 


(J  >/a(s)b(s)  ds)  - [yQ  /^+  xoeo'slnh^{  /a( 


s)b(s)  ds) } 


+ Ml  1 ds{a("s')b(s)  d^alsj]  }sinh(|  ^(ajb^)  da)  x(s)ds.  (U) 

0 


The  value  of  the  Volterra  integral  equation  (30),  however,  is  not  so  much  for 

direct  computation  as  it  is  for  suggesting  an  approximation,  the  so-called  Liouville- 

r 32 1 

Green  approximation  (see  OLVER  ),  to  the  solution  of  the  X force-level  equation 
(11).  If  the  appropriate  fractional  power  of  the  relative  effectiveness  is  "slowly 
varying,"  then  by  (21)  we  would  expect  that  |F(t)|  <<  1 so  that  we  could  drop  the 
integral  term  in  (30)  to  obtain 


X(t)  - xQ  cosh  (t-tq)  - [y0/a0/b0+  x0e0]sinh(T-T0) » 


(32) 


where  X(t)  denotes  the  Liouville-Green  approximation.  A theoretical  error  analysis 
[321  [331 

appears  in  Olver  (see  also  OLVER  ),  although  we  will  not  pursue  this  matter 
further  here.  We  observe  that  F(t)  a 0 V t i Tq  implies  that  while  X(t)  i 0 we 
have  X(t)  i X(t)  and  similarly  when  F(t)  & 0.  As  we  shall  see  below,  such  cases 
in  which  F(t)  is  always  i 0 or  * 0 are  readily  encountered  in  applications.  In 
terms  of  the  original  variables  t and  x,  we  have 


x(t)  * 


a(t)/a. 


-|1 


b(t)/br 


{Xq  cosh 


(j  /a(s)b(s)  ds)  - [yQ  + x0e0]sinh(j  /a(s)b(s)  ds)},  (33) 
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which  we  may  also  write  in  terms  of  the  average  combat  intensity. 


r 


(l/t)  /a(s)b(s)  ds , as 


x(t)  - 


a(t)/a 


-il 


b(t)/br 


*{xQ  cosh(/a(t)b(t)  t)  - [yQ/  x0e0Jsinh(/a(t)b(t)  t)}. 


(34) 


6.  Some  New  Victory-Prediction  Conditions. 

In  this  section  we  show  that  much  valuable  information  about  the  course  of 
combat  (for  example,  force-annihilation  prediction)  may  be  obtained  directly  from 
Liouville's  normal  form  (20)  without  making  any  kind  of  approximation.  Let  us  assume 
that  the  model  (7)  holds  for  all  time  and  consider  a fixed-force-ratio-breakpoint 
battle  (see  Section  2).  In  order  to  insure  that  the  battle  terminates  in  finite  time, 
we  must  make  certain  technical  assumptions  about  the  mathematical  nature  of  the  attri- 
tion-rate coefficients  in  (7).  For  the  reader's  convenience,  we  list  in  Table  I the 
principal  such  conditions  that  we  use  in  subsequent  developments.  We  observe  that  all 
these  conditions  hold  for  constant  attrition-rate  coefficients.  Let  us  also  observe 
that  Conditions  (B)  and  (ND)  imply  Condition  (C). 

We  record  here  for  future  reference  Taylor  and  Parry 's^^'  generalization  of 
Theorem  1 to  cases  of  variable  attrition-rate  coefficients  (see  Note  14). 


THEOREM  2:  Assume  that  Conditions  (B)  and  (ND)  hold.  Then  x^/y^  < ^aQ^Q 
implies  that  the  X force  will  lose  a fixed-force-ratio-breakpoint  battle 
in  finite  time. 


PROOF:  Introducing  the  force  ratio,  u = x/y,  we  have 

du/dt  - b(t)u2-a(t)  with  u(t=0)  * uQ  = xQ/yQ.  (35) 

Let  u+(t)  » /a(t)/b(t)  - /R(t)  denote  the  positive  root  of  the  quadratic  equation 
b(t)u2-a(t)  * 0,  and  observe  that  du/dt  < 0 for  any  u < u+(t)  (see  Figure  2 of 
Taylor  and  Parry^33^).  Condition  (ND)  means  that  u+(t)  is  nondecreasing.  It  is 
readily  shown  that  du/dt (t^O)  < 0 and  u+(t)  nondecreasing  imply  that  4u/dt(t)  < 0 
for  all  t i 0 (see  pp.  526-527  of  Taylor  and  Parry ^^).  Consequently,  from  (35)  we 
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TABLE  I.  Assumptions  about  Attrition-Rate  Coefficients. 


t t 

CONDITION  (A):  j a(s)ds  and  [ b(s)ds  are  bounded 


for  all  finite  t i tQ. 


i 

i 


CONDITION  (B):  lim  | b(s)ds  - 4~. 


c 

CONDITION  (C):  lim  /a(s)b(s)  ds  - 4«. 
t-*-H»  J 


CONDITION  (ND) : R(t)  is  nondecreasing  on  [0,+"). 
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r 


see  that  xQ/y0  * •'ag/bg  and  Condition  (ND)  yield  that  du/dt(t)  < 0 V t S 0.  It 
then  remains  to  be  shown  that  X's  breakpoint  force  ratio  is  reached  in  finite  time 
(see  Note  15).  Let  us  observe  that  R(t)  nondecreasing  and  0 < a(t),b(t)  < +*V  t (0,+“>) 
• aQ  < -M»  and  bQ  > 0.  From  Condition  (ND)  and  the  fact  that  du/dt(t)  < 0V  t a 0, 


we  have 


Hence 


£(t)  - b(t){'J2-(t(t»  * ^b0u2-a0)  - ^ ^(t-0). 


“(t) ' °o+  \ (dr)dt  ‘ ',o-fb7dir<t'0>  I b(*)d,> 


so  that  lim  b(s)ds  ■ +*  implies  that  u(t)  goes  to  il  i 0 in  finite  time 

t—  JQ 

We  observe  that  Theorem  2 says  that  the  X force  will  be  annihilated  in  finite 


Q.E.D. 


time  (i.e.  u(t«t  ) 
a 


X f 

0 for  tfl  finite  when  Uy  ■ 0). 


We  will  now  deduce  some  new  results  that  many  times  allow  us  to  predict  such  a 
battle's  outcome  from  the  initial  conditions.  Let  U = -(l/X)dX/dr  so  that  (20) 
becomes  the  Riccati  equation 
dU 


dr  ” U2  - 1 - F(t) 


with 


y /a 

U(T-Tn)  - — Ar+En.  (36) 
0 x0/  b0  0 


We  observe  that 


U(t)  - Z R1/2  + e(t). 


If  F(t)  « 0 V t a tq,  then  dU/dt  * U2-l,  and  U(t-tq)  - UQ  > 1 


(37) 


lim  U(t)  “ +<» 


for  t finite.  We  observe  from  (22)  and  the  definition  of  tn  that  R(t(t))  and 

a 

e(t)  are  finite  for  all  t a Tq  > 0 so  that  by  (37)  the  X force  must  be  annihilated 
in  finite  time,  since  by  Conditions  (A)  and  (C)  i(t)  is  a strictly  increasing  func- 
tion V t e [tQ,+")  with  range  [0,+®).  If  we  additionally  assume  that  dR/dt  *•  0, 

then  du/dt(t)  < 0 for  all  t e [0,t*]  with  lim  u(t)  » 0 so  that  the  X force  will 

t-*tZ 
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lose  any  f lxed-force-ratio-breakpoint  battle  In  finite  time.  The  proof  is  by 
contradiction:  dR/dt  * 0 ■»  u+(t)  nonincreasing  =>  i f there  exists  t.,  such  that 
dx/dt(t^)  i 0,  then  du/dt(t)  i 0 V t i t^  ° impossible  to  have  u -*■  0.  Thus,  we 
have  proved 


THEOREM  3:  Assume  that  Conditions  (A)  and  (C)  hold.  If  F(t)  S OVt  i t^, 

then  — (l-e„) < / -r^-  implies  that  the  X force  will  be  annihilated  in  finite 
y0  0 /b0 

time.  Furthermore,  if  dR/dt  < OVt  4 then  the  X force  will  lose  any 


fixed-force-ratio-breakpoint  battle  in  finite  time. 


If  F(x)  i 0 and  dR/dt  iOVtir,  then  dU/dr  S U2  - 1 and  0,  < 1 « 3 t 

U <j 

finite  such  that  U(t3)  = 0 with  U(t)  > 0 V t e [Tq,^).  Since  dR/dt  i 0 =»  e(t)  i 0, 
we  see  from  (37)  that  Y must  be  annihilated  in  finite  time,  since  f(t)  is  a one- 
to-one  mapping  of  [tp,+“)  onto  [0,+“>).  By  observing  that 


we  see  that  du/dt(t)  > OV  t i 0 as  long  as  u > 0.  Hence  u(t)  is  strictly  increasing 

cl  SL 

for  te  [O.ty)  and  lim  u(t)  = +“>  for  t^  finite  so  that  Y will  lose  any  battle 
t>t- 

f 3 

with  ux  > Uq.  Thus,  we  have  proved 


THEOREM  4:  Assume  that  Conditions  (A),  (B) , and  (ND)  hold.  If  F(t)  i 0V  Ti  t^, 


then  — (1— e ) > 
y0 

breakpoint  battle  in  finite  time. 


/5 

/bo 


implies  that  the  Y force  will  lose  a f ixed-force-ratio- 


By  considering  Lioville's  normal  form  (24)  for  the  Y force-level  equation, 
one  may  similarly  prove  Theorems  5 and  6. 

THEOREM  5:  Assume  that  Conditions  (A)  and  (C)  hold.  If  G(t)  £ OVt^Tq, 
then  — > Pp  (l+ert)  implies  that  the  Y force  will  be  annihilated  in  finite 

y0  Vb0  0 

time.  Furthermore,  if  dR/dt  i 0 V x i (i.e.  Condition  (ND)  holds),  then  Y 
will  lose  any  fixed-force-ratio-breakpoint  battle  in  finite  time. 


Cfl  >-■ 


THEOREM  6:  Assume  that  lim 


t-*-H»  J 0 


II 


a(s)ds  = +*»  and  that  Condition  (A)  holds.  If 


X /EL 

G(t)  i 0 and  dR/dt  < OV  T i t.,  then  — < Jrp-  (1+e-)  implies  that  the  X 

0 yo  y bo  ° 

force  will  lose  a fixed-force-ratio-breakpoint  battle  in  finite  time. 


Although  F(t)  and  G(t)  are  lated  by  (29),  for  all  the  attrition-rate 
coefficients  that  we  have  so  far  considered,  F(t)  i 0 *>  G(t)  £ 0.  Also,  in  all  these 
cases  F(t)  i 0 *>  dR/dt  i 0.  Since  xQ/y0  < '^Tb^  (l+eQ)  **  (xc/y0^1-G0^  < 
and  (x^/y^) (1-eq)  > ^aQ/b0  =»  xQ/y0  > /a0/b~  (l+eQ) . we  see  that  Theorem  6 is  stronger 
than  Theorem  3 and  Theorem  4 stronger  than  Theorem  5.  Under  these  conditions.  Theorems 
4 and  6 provide  no  additional  information  on  battle  outcome  over  that  contained  in 
Theorems  3 and  5. 


7.  Application  to  General  Power  Attrition-Rate  Coefficients. 

We  now  apply  the  above  results  to  combat  between  two  homogeneous  forces  modelled 
by  (7)  with  the  general  power  attrition-rate  coefficients  (8).  We  distinguish  between 
two  cases:  (I)  power  attrition-rate  coefficients  with  no  offset  (i.e.  A=0),  and 

(II)  power  attrition-rate  coefficients  with  the  same  parity  (i.e.  p = v)  and  positive 
offset.  In  order  to  invoke  Theorems  3 through  6,  we  must  have  C > 0 for  the  general 
power  attrition-rate  coefficients  (8)  (cf.  (19)).  In  order  that  Condition  (A)  be 
satisfied,  we  must  have  p,v  > -1;  and  then  Conditions  (B)  and  (C)  hold  (see  Table  I). 

7.1.  Power  Attrition-Rate  Coefficients  with  No  Offset. 

In  this  case  we  have 

a(t)  - ka(t+C)W  and  b(t)  - k^t+C)'',  (38) 

so  that 

dR/dt  i 0 «»  v i v.  (39) 

Theorem  2 then  yields 

COROLLARY  2.1:  For  the  power  attrition-rate  coefficients  (38)  with  p i v and 
C > 0,  — < J implies  that  the  Y force  will  win  a f ixed-force-ratio- 

y0  Vbo 
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breakpoint  battle  in  finite  time.  In  particular,  the  X force  will  be 

annihilated  in  finite  time  when  u^  = 0. 

In  preparation  for  invoking  Theorems  3 through  6,  we  compute 

Fm  = (P-v)  (3u+v+4) 

F(  ' 4(p+v+2)ztz  ’ 


and 


since  (13)  with  (38)  yields 


. . (v-p) (p+3v+4) 

G(T)  “ 4(u+v+T)/t/  * 


T-T(t)  =^#(t+C)(^2)/2 


(p+v+2) 


,(p+v+2)/2 


with  Tq  » ( 2 *iT^k^7  ( p+v+2 ) )C'H  ' v ' > 0.  We  observe  that 

F(t)  i 0 and  G(t)  S 0V  i > 0 «»  p i v. 


(AO) 


(Al) 


(A2) 


(A3) 


Thus,  as  discussed  above  at  the  end  of  Section  6,  Theorems  A and  6 provide  no  additional 
information  on  force  annihilation  over  that  contained  in  Theorems  3 and  5.  Furthermore, 
Theorem  3 is  stronger  than  Theorem  2,  since  e^,  shich  is  given  by  (22)  evaluated  at 
t * 0,  is  i 0.  Hence,  we  omit  any  corollary  to  Theorem  3.  As  a corollary  to 
Theorem  5 we  have 


COROLLARY  5.1:  For  the  power  attrition-rate  coefficients  (38)  with  p ^ v 
and 


C > 0,  — > U+  (y~v?  c-(w+v+2)/2}  impiies  that  the  X force 

y0  b0  A/k  k 


S 

will  win  a fixed-force-ratio-breakpoint  battle  in  finite  time.  In  particular, 


thd  Y fotce  will  be  annihilated  in  finite  time  when  uv 


+“>. 


Let  us  also  write  the  X-victory-prediction  condition  of  Corollary  5.1  as 


- f(c). 

y„  V k,  Ak^ 


(AA) 


Observing  that  lim  f(C)  “ +®  for  p > v,  we  see  that  the  victory-prediction  condi- 
C-K)+ 

tlons  of  Corollaries  2.1  and  5.1  become  stronger  as  C decreases  and  are  meaningless 
for  C » 0. 
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We  now  show  that  other  considerations,  however,  readily  yield  force-annihilation 
conditions  when  C = 0.  [Unfortunately,  without  further  analysis  these  results  are 
limited  to  force-annihilation  prediction  and  do  not  apply  directly  to  fixed-force-ratio- 
breakpoint  battles  with  u^  < +°°  and  u^  > 0.]  The  solution  to  (7)  with  the  power 
attrition-rate  coefficients  (38)  may  be  written  for  C = 0 as  (see  Note  16) 

x(t)  = (pP/(2p))(Cx+CY}AB(T)+  (pP/(2^p)){Cx-Cy}BB(T),  (45) 


where  AD  and  B denote  the  generalized  Airy  functions  of  the  first  and  second  kinds 
of  order  8 (see  SWANSON  and  HEADLEY ) , p = (u+1) /(p+v+2) , 8 = (v-p)/(y+l), 

Cx  = x0  r(l-p) , = y0^7^b  F(p)  (viTk^/(u+\H-2))1_2p,  and  T = (/klc^/(p+l))2ptli+1. 

Observing  that  A (O,  B (O  > 0V(  i 0,  lim  A (?)  = 0,  and  lim  B (£)  = -H», 

V ^->+0O  V 

we  see  from  (45)  that  lim  x(t)  = if  and  only  if  the  coefficient  of  B (T)  is 

t-*~H»  6 

negative  (i.e.  Cx  < C^) . Hence,  we  may  conclude  (cf.  the  development  of  (3)  from  (5)) 


PROPOSITION  1:  For  a f ight-to-the-f inish  modelled  with  the  power  attrition- 
rate  coefficients  (38)  and  C = 0 , the  X force  will  be  annihilated  in  finite 
time  if  and  only  if  ? • 

For  C > 0,  (45)  does  not  take  nearly  such  a convenient  form  (see  Taylor  and  Brown  ^°^), 

and  the  analogue  of  Proposition  1 for  C > 0 would  involve  the  generalized  Airy  func- 
tions and  their  derivatives  evaluated  at  a positive  argument.  Thus,  the  prediction  of 
force  annihilation  by  this  approach  for  C > 0 would  require  tabulations  of  higher 
transcendental  functions.  Such  tabulations  do  not  currently  exist  (see  reference  50). 
Thus,  we  see  the  usefulness  of  the  "strong"  annihilation  conditions  given  by  Corollaries 
2.1  and  5.1,  which  contain  only  elementary  functions:  "exact"  annihilation-prediction 
conditions  involve  higher  transcendental  functions. 

7.2.  Power  Attrition-Rate  Coefficients  with  the  Same  Parity  and  Positive  Offset. 

In  this  case  we  have 


a(t)  - ka(t+C)‘ 


b(t)  = kb(t4C+A)W, 


with  A,C  > 0.  It  follows  that 


1 


dR/dt  i 0 «»  v i-  0. 


(47) 


Then  Theorem  2 yields 

COROLLARY  2.2:  For  the  offset  power  attrition-rate  coefficients  (46)  with 
y i 0 and  A,C  > 0,  — < J r~  implies  that  the  Y force  will  win  a fixed- 

yo  y bo 

force-ratio-breakpoint  battle  in  finite  time. 

Let  us  now  see  what  information  is  yielded  by  Theorems  3 through  6.  It  is  more 
convenient  to  express  F and  G in  terms  of  the  old  time  variable  t (see  (27)  and 
(28)).  Then 


16^1^  ( t4C ) V+2  ( t+C+A)  W+2 


G(x)  = 


16k  k ( t+C  ) *J+2  ( t+C+A)  y+2 
a b 


We  observe  that 


F(t)  i 0 and  G(t)  sOVt  > 0 « p £ 0,  (50) 

so  that  Theorems  4 and  6 again  provide  no  additional  information  over  that  contained 
in  Theorems  3 and  5.  Furthermore,  we  may  omit  consideration  of  Theorem  3,  since  it  is 
implied  by  Theorem  2.  As  a corollary  to  Theorem  5 we  have 

COROLLARY  5.2:  For  the  offset  power  attrition-rate  coefficients  (46)  with 


U S 0 and  A.C  > 0,  — > U + 

y~  v b_ 


4/aQb0  C (C+A) 


} implies  that  the  X 


force  will  win  a fixed-force-ratio-breakpoint  battle  in  finite  time. 
Let  us  also  write  the  X-victory-prediction  condition  of  Corollary  5.2  as 


<i«/c)-“,2tjE-o^5ST-  g(C). 


We  observe  that  lim  g(C)  = +00,  and  we  again  see  that  the  victory-prediction  conditions 
C-K)+ 

of  Corollaries  2.2  and  5.2  become  stronger  as  C decreases  and  become  meaningless  for 
C = 0. 


I 


8.  Discussion. 
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In  his  classic  1914  paper  1 , Lanchester  assumed  that  the  combatants'  fire 

effectivenesses  (as  expressed  by  the  Lanchester  attrition-rate  coefficients)  were 
constant  over  time  and  deduced  his  famous  square  law  (2) , which  allows  one  to  tradeoff 
quality  versus  quantity  of  weapon  systems  by  means  of  the  condition  for  equality  of 
"fighting  strengths" 

xQ/y0  = *4/b,  (52) 

where  a and  b denote  constant  attrition-rate  coefficients  (see  Section  2).  Thus, 
we  see  that  equality  of  Lanchester-type  fighting  strengths  (see  Note  17)  depends  on 
two  parameters:  the  initial  force  ratio  and  the  relative  effectiveness.  When  the 
timing  of  military  actions  is  considered,  we  add  a third  parameter,  the  intensity  of 
combat,  to  this  list  of  significant  combat  parameters.  In  the  paper  at  hand,  we 
extended  these  well-known  constant-coefficient  results  to  battles  between  two  homogeneous 
forces  with  temporal  variations  in  the  fire  effectivenesses. 

No  such  simple  relationship  like  the  square  law  (2) , which  yielded  (52) , holds 
in  general  for  variable  attrition-rate  coefficients  (see  Note  18) . By  transforming 
the  independent  variable  to  normalize  the  battle's  time  scale  by  the  intensity  of  com- 
bat, we  found  that  the  course  of  combat  depends  on  two  weapon  system  parameters:  (I) 
relative  fire  effectiveness,  R(t)  = a(t)/b(t),  and  (II)  intensity  of  combat,  I(t)  = 
>/,a(t)b(t) . Moreover,  when  the  temporal  variations  in  relative  fire  effectiveness 
R(t)  follow  a regular  pattern  (e.g.  R(t)  nondecreasing),  the  battle's  outcome  can 
many  times  be  predicted  from  the  battle's  initial  conditions.  To  obtain  such  results, 
we  considered  Liouville's  normal  form  for  the,  for  example,  X force-level  equation 
and  found  that  it  not  only  yields  new  battle-outcome-prediction  conditions  but  also 
suggests  an  approximation  to  the  time  history  of  the  X force  level  x(t).  As  seen 
from  (21),  Liouville's  normal  form  (20)  introduces  second  order  conditions  (e.g.  the 
second  derivative  of  relative  fire  effectiveness  with  respect  to  transformed  time)  into 
Lanchester-type  combat  analysis. 


i: 
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The  new  victory-prediction  results  (see  Theorems  3 through  6)  that  we  have 
developed  here  are  complementary  to  those  of  Taylor  and  Parry (see  Theorem  2).  This 
complementary  nature  is  seen  by  observing  that  under  the  appropriate  comparable  condi- 
tions we  have 


from  Taylor  and  Parry's  Theorem  2 (4^- i 0)  : 

dt 

from  our  new  results.  Theorem  5 (G(t)  £ 0) : 


Y will  win  if 

X will  win  if 


where 


(53) 


For  both  special  cases  of  general  power  attrition-rate  coefficients  considered  in 
Section  7 above  we  had  dR/dt  S 0 « F(t)  k 0 G(t)  t 0 e i 0,  Although  these  if- 
and-only-if  statements  do  not  hold  in  general,  they  do  hold  for  these  particular  coeffi- 
cients. In  both  cases,  we  observe  that  for 


/aot!on 

b0  *0  /bo  0 


(54) 


we  cannot  say  by  this  approach  who  will  be  the  loser  of  the  f ixed-force-ratio-break- 
point  battle  (see  Figure  3) . We  observe  that  force  annihilation  is  a special  case  of 
these  outcome-prediction  conditions.  From  both  (54)  and  Figure  3,  we  see  that  there 
is  a "gap"  in  these  victory-prediction  conditions  (i.e.  Theorems  2 through  6).  The 
price  of  removing  this  "gap,"  however,  is  the  introduction  of  higher  transcendental 
functions  (see  Section  7 above  and  Taylor  and  Brown ^^).  Furthermore,  "exact"  results 
with  no  such  gap  in  the  victory-prediction  conditions  are  apparently  only  possible 
for  a fight-to-the-finish  in  which  one  side  or  the  other  is  to  be  annihilated. 

We  also  refined  a victory-prediction  result  of  Taylor  and  Parry (as  applied 
to  the  model  (7)  under  study)  by  adding  a restriction  on  the  attrition-rate-coefficients, 
i.e.  our  Condition  (B),  to  the  assumptions  of  Theorem  2.  Taylor  and  Parry  did  not 
note  that  such  a condition  must  be  assumed  in  order  to  insure  that  the  battle  will 
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Figure  3.  "Gap"  in  the  Victory-Prediction  Conditions  for  dR/dt  i 0,  F(t) 


terminate.  All  our  new  results,  i.e.  Theorems  3 through  6,  contain  some  such  restriction 
on  the  attrition-rate  coefficients. 

The  results  of  this  paper  may  be  used  in  parametric  analyses  (see  BONDER^)  of 

the  dynamic  combat  interactions  between  two  homogeneous  forces  with  time-(or  range-) 

dependent  weapon  system  capabilities.  Such  models  are  of  particular  interest  in  light 

of  the  work  of  S.  BONDER^’^  and  others  on  the  prediction  of  Lanchester 

attrition-rate  coefficients  from  weapon  system  performance  data  and  the  work  of 

[12] 

G.  CLARK  on  the  estimation  of  such  (time-dependent)  coefficients  from  Monte  Carlo 
simulation  output.  A further  discussion  of  applications  is  to  found  in  references  8 
and  43.  As  is  always  the  case,  however,  the  insights  gained  into  combat  dynamics  from 
such  Lanchester-type  models  are  no  more  valid  than  the  models  themselves. 

9.  Summary. 

In  this  paper  we  have  developed  insights  into  the  dynamics  of  Lanchester-type 
combat  between  two  homogeneous  forces  with  temporal  variations  in  weapon  system 
effectivenesses  by  considering  Liouville's  normal  form  for  the  X and  Y force-level 
equations.  Our  principal  results  were  some  new  outcome-prediction  results  that  com- 
plemented those  of  Taylor  and  Parry We  also  saw  that  Taylor  and  Parry's  victory- 
prediction  result  (as  applied  to  our  model  (7)),  Theorem  2,  must  be  refined  by  making 
an  additional  assumption  about  the  attrition-rate  coefficients.  Liouville's  normal 
form  also  suggests  an  approximation,  the  Liouville-Green-Lanchester  approximation,  to 
the  force-level  trajectories.  Additionally,  by  transforming  the  battle's  time  scale, 
we  saw  that  the  relative  fire  effectiveness  of  the  two  combatants  and  the  intensity 
of  combat  were  two  key  parameters  affecting  the  course  of  battle.  Our  new  victory- 
prediction  conditions  for  fixed-force-ratio-breakpoint  battles  were  applied  to  two 
special  cases  of  combat  modelled  with  general  power  attrition-rate  coefficients:  (I) 
power  attrition-rate  coefficients  with  no  "offset"  (modelling,  for  example,  two  weapon 
systems  with  the  same  miximum  effective  range),  and  (II)  offset  power  attrition-rate 
coefficients  with  the  same  parity  (modelling,  for  example,  weapon  systems  with  different 


maximum  effective  ranges  but  the  same  type  of  range  dependence).  We  saw  that  there 
was  a "gap"  in  the  range  of  initial  force  ratios  for  which  we  could  predict  the  outcome 
of  such  a fixed-force-ratio-breakpoint  battle.  The  price  that  one  has  to  pay  to 
remove  this  "gap"  was  discussed. 


NOTES 


1.  F.  W.  Lanchester  (1868-1946)  was  an  English  automotive  and  aeronautical  engineer. 

For  a brief  sketch  of  his  many  scientific  and  engineering  contributions,  see  McCLOSKEY 
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In  acknowledgment  of  his  contribution  to  operations  research  (again,  see 
reference  29)  the  Operations  Research  Society  of  American  annually  awards  the  Lanchester 
Prize  (see  p.  113  of  reference  34)  "for  the  paper  on  operations  research  judged  to  be 
the  best  of  the  calendar  year." 


2.  See,  for  example,  Bonder  and  Farrell^81,  Bonder  and  Honig^91,  Brackney  1 33  1 , 
DEITCHMAI^141,  MORSE  and  KIMBALL^ 30 ^ , SCHAFFER^391,  Taylor  and  Parryt531,  WALLIS^521, 


and  Weiss 


[54,55] 


3.  As  work  by  Bonder  and  Farrell^81,  Taylor^431,  and  Taylor  and  Brown^3 4^1  shows,  the 
infinite-series  solution  to  variable-coefficient  equations  by  itself  provides  little 
information  about  battle  outcome  because  of  its  complexity. 

4.  In  his  well-known  survey  paper  on  the  Lanchester  theory  of  combat,  Dolansky^331 
suggested  the  development  of  outcome  predicting  relations  without  solving  in  detail 
and/or  computing  force-level  trajectories  as  one  of  several  problems  for  future  research 
The  work  at  hand  is  a step  towards  this  problem's  resolution  (see  also  Taylor  and 
Parry1331  and  Taylor1481)- 
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5.  Scientific  verification  of  Lanchester-type  models  (as  with  any  combat  model)  is 
still  an  unresolved  question  (see  Bonder ^).  Although  there  have  been  numerous  attempts 
to  compare  the  theoretical  implications  of  such  models  (invariable  quite  simple,  con- 
stant-coefficient ones)  with  empirical  (i.e.  historical)  evidence  (see , for  example, 
ENGEL^16^,  HELMBOLD^19,20^ , SCHMIEMAN*4^ , WEISS ^ 54 ’ 57 ^ , and  WILLARD[58^) , the  results 
have,  unfortunately,  been  inconclusive,  with  far  from  universal  agreement  as  to  their 
correct  interpretation.  The  historical  data  base  is  apparently  not  rich  enough  in 

detail  to  permit  a definitive  answer  to  the  scientific  validity  of  Lanchester-type 

[21] 

models  (see  HELMBOLD  ),  since  nations  fight  wars  for  other  reasons  than  to  collect 
combat  data. 

6.  H.  K.  WEISS^8^  has  pointed  out  that  Lanchester,  an  Englishman,  was  anticipated  (in 
qualitative  but  not  quantitative  terms)  in  1905  by  Bradley  A.  Fiske  (then  Commander 
but  later  Rear  Admiral,  USN) , an  American.  For  a sketch  of  the  life  and  accomplish- 
ments of  Bradley  Allen  Fiske  (1854-1942) , see  pp.  298-299  of  reference  31.  J.  ENGEL^7^ 
subsequently  showed  that  Fiske's  verbal  model  is  equivalent  to  a system  of  difference 
equations  (in  contrast  to  Lanchester's  differential  equations)  and  examined  some  of 

the  mathematical  consequences  of  these  Fiske-type  equations  of  warfare. 

7.  The  equations  (1)  are  only  valid  for  x,y  > 0.  The  first,  for  example,  becomes 
dx/dt  * 0 for  x - 0.  Moreover,  there  is  far  from  universal  agreement  as  to  which 
variables  are  significant  and  can  be  used  to  predict  the  outcome  of  the  combat  process 
For  some  other  views,  see  HAYWARD^ 18 ^ and  LIDDELL  HART[27^. 

8.  The  influential  19th-century  German  military  philosopher,  Carl  von  Clausewitz 
(1780-1831),  stated  in  his  classic  work  On  War  (Vom  Kriege)  (see  p.  276  of  reference 
13),  "The  best  Strategy  is  always  to  be  very  strong,  first  generally  then  at  the  deci- 
sive point.  ...There  is  no  more  imperative  and  no  simpler  law  for  Strategy  than  to 
keep  the  forces  concentrated." 
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9.  As  pointed  out  by  Taylor  and  Parry^"*^,  the  entire  subject  of  modelling  battle 

termination  is  a problem  area  in  contemporary  defense  planning  studies.  There  is  far 

[47] 

from  universal  agreement  on  this  topic  (see  TAYLOR  for  further  references). 

10.  Before  the  mid-1960's  the  use  of  Lanchester-type  models  in  defense  planning  studies 

[41 

was  hampered  by  the  inability  to  predict  the  attrition-rate  coefficients  (see  Bonder  ). 
Thus,  two  significant  accomplishments  for  the  Lanchester  theory  of  combat  in  the  1960's 
were  (1)  the  development  of  methodology  for  the  prediction  of  Lanchester  attrition-rate 
coefficients  from  weapon  system  performance  data  by  S.  Bonder  and  others  (see 

BARFOOT^,  Bonder  and  Farrell^,  and  KIMBLETON  ) , and  (II)  G.  Clark's^12^  develop- 
ment of  methodology  for  the  (maximum  likelihood)  estimation  of  such  coefficients  from 
Monte  Carlo  simulation  output.  Both  these  developments  and  others  (see  references 

10.  37,  and  38)  have  facilitated  the  application  in  defense  planning  studies  of  models 
such  as  (7)  and  its  generalization  to  combat  between  heterogeneous  forces  (see  refer- 
ence 8) . 

11.  The  modelling  roles  of  A and  C are  discussed  in  Taylor  and  Brown . 

12.  Taylor  and  Brown^"*^  give  an  example  of  nonexistence  of  a solution  to  (7)  for 

the  power  attrition-rate  coefficients  (8)  with  A = 0 (see  Section  4 of  reference  50). 

13.  The  term  quasi-autonomous  was  coined  by  Taylor^^  (see  also  TAYLOR^^)  to  denote  a 

system  of  differential  equations  transformable  to  an  autonomous  system  (see , for 

r 35 1 

example,  p.  163  of  PETROVSKI  ) by  a change  of  the  time  scale.  Special  cases  of 

[8] 

such  Lanchester-type  equations  have  been  considered  by,  for  example,  Farrell1  J and 
[42] 

TAYLOR  . More  general  (possibly  nonlinear)  quasi-autonomous  Lanchester-type 
equations  have  been  studied  by  Taylor (see  also  Note  4 of  Taylor  and  Brown ^^). 
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1A.  A similar  theorem  for  a more  general  model  with  supporting  fires  was  given  by 
Taylor  and  Parry , but  they  did  not  observe  in  reference  51  that  certain  additional 
conditions  must  be  assumed  to  insure  that  the  battle  will  terminate.  In  other  words, 
an  analogous  result  originally  given  by  Taylor  and  Parry for  such  a linear,  vari- 
able-coefficient model  with  supporting  fires  is  not  true  in  general  without  certain 
restrictions  on  attrition-rate  coefficients  being  added. 

15.  This  point  was  overlooked  by  Taylor  and  Parry^^  (see  Note  14  above). 

16.  The  substitution  s = K/k^/k^  j a(o)do  where  K = (/^k^/ (y+1) ) 2p  1 transforms 
the  X force-level  equation  (11)  with  power  attrition-rate  coefficients  (38)  and 

Q 

C = 0 into  d2x/dx2-sx  = 0,  which  is  readily  recognized  as  the  generalized  Airy 

[All 

equation  (see  Swanson  and  Headley  ). 

17.  The  determination  of  equality  of  fighting  strengths  will  be  affected  by  the 
battle  termination  model  used  in  the  operational  definition  of  the  concept  of  fighting 
strength.  For  a discussion  of  related  matters,  see  reference  18  in  which  P.  Hayward 
examines  the  factors  that  determine  combat  effectiveness. 

18.  As  apparently  first  observed  by  B.  0.  Koopman^^,  a "square-law"  relationship 
holds  for  the  quasi-autonomous  case  in  which  a(t)/b(t)  = constant  (see  also  refer- 
ence 50) . 
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